Single particle relaxation time versus transport scattering time in a 2D graphene layer 
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We theoretically calculate and compare the single-particle relaxation time (r s ) defining quantum 
level broadening and the transport scattering time (r t ) defining Drude conductivity in 2D graphene 
layers in the presence of screened charged impurities scattering and short-range defect scattering. 
We find that the ratio Tt/r s increases strongly with increasing fc^Zi and k where kw, Zi, and n 
are respectively the Fermi wave vector, the separation of the substrate charged impurities from the 
graphene layer, and the background lattice dielectric constant. A critical quantitative comparison of 
the Tt/r 3 results for graphene with the corresponding modulation-doped semiconductor structures 
is provided, showing significant differences between these two 2D carrier systems. 

PACS numbers: 81.05.Uw; 72.10.-d, 72.15.Lh, 72.20.Dp 



I. INTRODUCTION AND BACKGROUND 

In an impure disordered conductor (either a metal or a 
doped semiconductor) scattering by static impurities and 
defects, in general, leads to two distinct momentum relax- 
ation times; the scattering life time or the transport re- 
laxation time (denoted by r t in this paper, but often just 
called r) and the quantum lifetime or the single-particle 
relaxation time (denoted by r s in this paper, but often 
also called r q ). Although r t and r s both arise from impu- 
rity scattering in the metallic regime, they are in general 
distinct and unique with no direct analytical relationship 
connecting them, except for the simple (and often unre- 
alistic) model of completely isotropic s-wave zero-range 
impurity scattering when they are equal. In particular, 
Tt determines the conductivity, a — ne/x oc Tt where n is 
the carrier density and [i the mobility, whereas r s deter- 
mines the quantum level broadening, T = H/2t s , of the 
momentum eigenstates (i.e. the band states defined by 
a momentum fik with k as the conserved wave vector). 
The difference between r t and t s arises from the subtle 
effect of the wave vector dependent impurity potential 
Ui (k) which distinguishes between momentum scattering 
in ail directions contributing to r.T 1 and transport relax- 
ation TfT 1 which is unaffected by forward scattering (or 
more generally, small-angle scattering)^. 

In this paper we study theoretically r t and r s in the 2D 
graphene layer due to long-range and short-range impu- 
rity scattering, finding interesting behavior in the ratio 
r t /r s as a function of impurity location, carrier density, 
and the system environment (e.g. the background lat- 
tice dielectric constant of the substrate). We compare 
graphene r t /r s with the corresponding situation in the 
extensively studied 2D semiconductor systems, finding 
significant difference. 

In 3D metals or semiconductors, the effective impurity- 
electron interaction potential is almost always short- 
ranged. This is true even when the bare impurity poten- 
tial is long-range Coulombic since carrier-induced elec- 
tronic screening of the charged impurities is highly effec- 
tive. This strongly screened short-range impurity poten- 



tial in 3D systems (as well as in 2D systems when elec- 
trons and impurities are not spatially separated) has led 
to the almost universal theoretical adoption of the uncor- 
related white-noise zero-range impurity potential model 
for studying the impurity scattering effect on "metallic" 
transport properties. Since such a white-noise disorder, 
by definition, leads to isotropic impurity scattering dom- 
inated entirely by large-angle scattering, r t = t s in or- 
dinary metals and semiconductors for such a zero-range 
disorder model as forward scattering simply plays no spe- 
cial role here. Graphene is, however, qualitatively differ- 
ent even for this short-range white noise disorder model 
due to its chiral sublattice symmetry which suppresses 
backward (i.e. a scattering induced wave vector change 
by 2Uf from +kp to —kp) scattering^, thus introducing 
an intrinsic chiral preference for forward scattering over 
backward scattering. As we show later in this paper, due 
to this suppression of backward scattering, the zero-range 
white-noise disorder potential model leads to r t = 2r s in 
graphene in contrast to the r t = t s in ordinary metals 
and semiconductors. This is due to the importance of 
kp (rather than 2k f) scattering in dominating graphene 
transport properties whereas the small-angle scattering 
always dominates t s . 

For long-range impurity scattering, however, it has 
been shown^ that r* can exceed r s by a large factor, 
and in high-mobility 2D modulation-doped semiconduc- 
tor structures, where the predominant scattering is the 
small-angle scattering due to charged impurities located 
far from the 2D electron layer, in general r t 3> t s . In 
this paper, we theoretically calculate r t /r s ( as well as 
Tt and t s individually) for 2D graphene electrons (and 
holes) due to both long-range and short-range impurity 
scattering, taking into account the dependence on the 
impurity location. We find that there are significant dif- 
ferences between 2D graphene (with its linear chiral dis- 
persion) and 2D semiconductors, showing that r t /r s in 
general could also be very large for graphene although 
the behavior is quite different from the corresponding 
2D semiconductors. 

Our work is restricted entirely to extrinsic or doped (or 
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gated) graphene where the Fermi level is away from the 
charge neutral Dirac point (taken to be the energy zero) . 
There has been extensive recent theoretical work 3 -^^^ 
on the charged impurity scattering limited transport 
properties in gated extrinsic graphene, but the single- 
particle quantum relaxation time r s (and its relation to 
the transport scattering time r t ) has not been discussed. 
We note that while r t has been experimentally studied re- 
cently in extrinsic graphene^ through measurements of 
density-dependent conductivity, the corresponding t s has 
not yet been directly experimentally measured. The most 
direct way of measuring t s experimentally is by measur- 
ing the quantum level broadening Y = h/2r s through the 
Dingle temperature To = Y jixk-B- The Dingle tempera- 
ture can be measured in low-field SdH magnetoresistance 
oscillations as has been done extensively in 2D semi- 
conductor structures^. Any measurement of the single- 
particle level broadening, e.g. magnetization, also gives 
an estimate for t s . Experimental comparisons of Tt / t s for 
various system parameters (e.g. carrier density, impurity 
location, materials parameters such as effective mass and 
lattice dielectric constant) have been carried out in semi- 
conductor based 2D systems, but not in graphene since 
measurements of t s are unavailable in graphene. Exten- 



sive measurements of Tt as a function of carrier density, 
i.e. gate voltage, however, have been carried oub^ in 
graphene recently. 

In our theoretical calculations (and in the numerical 
results) we have used two complementary models of bare 
disorder: Long-range random charged impurity centers 
at or near the graphene-substrate interface and short- 
range model white noise disorder in the graphene layer. 
We also include screening (by the graphene carriers them- 
selves) within random phase approximation (RPA) in the 
theory. In Sec. II we describe theoretical details and 
present our results. In Sec. Ill we discuss the relevance 
and significance of our theory in the context of transport 
experiments in gated or doped graphene, and provide a 
conclusion. 



II. THEORY AND RESULTS 

In the relaxation time approximation the graphene 
transport scattering time Tt by randomly distributed im- 
purity centers is given by 3 ^ 
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where rii is the concentration of the impurity center, 
q = |k — k'|, 6>kk' is the scattering angle between the 
scattering in- and out- wave vectors k and k', A, A' = ±1 
denote the band indices, e A k = A'yjlcj is a single par- 
ticle energy (7 = hv F being the graphene band veloc- 
ity), and F AA /(k, k') is the overlap of states and given by 
F\\i (k, k') = (l+AA'cos0 k)k ')/2. In Eq. Q V ei (q) is 
the matrix elements of the scattering potential between 
an electron and an impurity. For charged impurities, we 
use the Coulomb interaction V e i(q) = 2ire 2 /(nq), where 
k is the dielectric constant of surrounding materials, and 
for short-range point defect scatterers, V e %{q) = Vo, a con- 
stant. In Eq. [TJ e(q) is the static RPA dielectric (screen- 
ing) function appropriate for grapheneii, given by 



e(q) =l + « c (g)n(g), 



(2) 



where v c (q) = 2ire 2 /nq is the electron-electron Coulomb 
potential, and IL(q r ) is the polarizability function of 
graphene which is calculated to be 
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Since we consider elastic scattering we can neglect inter- 
band scattering processes (A' 7^ A). Then, we have the 



leading-order (in impurity disorder n») transport scatter- 
ing time 

L.^r^wf^.^ (4) 

T t 2-Kh Y J e(q) z 

where q = 2k F sin(0/2). 

From a many-body-theory viewpoint the single particle 
relaxation time t s can be calculated from the electron 
self energy of the coupled electron-impurity syste m 12 ' 13 . 
The electron self-energy due to the impurity scattering is 
given by 



E A (fc,u,) = ^ 



A'k' 



< l^(g) 2 | > 

e(q) 2 



F A v(k,k')G A (k», (5) 



where G A (k, cj) = l/(u> — e A k + i$) is the noninteraction 
Green's function. The single particle relaxation time is 
related to the imaginary part of the single particle self- 
energy function by 
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with the single particle quantum (impurity induced) level 
broadening Y s — lmS(k F , E F ), i.e. r s = h/2Y s . In the 
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leading order disorder approximation, we have the single 
particle relaxation time t s 



with 
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Thus, the only difference between the scattering time 
r t [Eq. and the single-particle relaxation time r s 

[Eq. ©] is the weighting factor 1 — cos 9 in the trans- 
port scattering time. The factor 1 — cos (9 weights the 
amount of backward scattering of the electron by the 
impurity. Small angle (or forward) scattering, where 
cos 9 w 1, is relatively unimportant in contributing to 
TjT 1 and contributes little to the resistivity. In normal 
2D systems the factor 1 — cos# obviously favors large 
angle scattering events, which are more important for 
the electrical resistivity. However, in graphene the large 
angle scattering is also suppressed due to the wave func- 
tion overlap factor 1 + cos 9. The transport scattering 
time thus gets weighted by an angular contribution fac- 
tor of (1 — cos 6){1 + cos 9), which suppresses both small- 
angle scattering and large-angle scattering contributions 
in the transport scattering rate. However, the single par- 
ticle relaxation time is weighted only by 1 -I- cos 9 term. 
Therefore, r t is insensitive to both small and large angle 
scatterings while r s is only sensitive to small angle scat- 
tering events. In fact, the dominant contribution to r t 
comes from cos 2 9 — 0, i.e. 9 = n/2 scattering, which is 
equivalent to |k— k'| = kp "right-angle" scattering in con- 
trast to the 2Uf back scattering in ordinary 2D systems. 
This leads to significant difference between graphene and 
parabolic 2D semiconductor systems with respect to the 
behavior of T{/t s . 

Using RPA screening function [Eq. j2])] at T = we 
calculate the scattering times due to charged impurities 
distributed completely randomly at the interface between 
graphene and the substrate with density rii C 
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In Eq. ([8|) Itc(x), Isc(x) are calculated to be 
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From Eqs. (fTTj) and (fT2)l we have 



for x < 1 

for x = 1 ■ (12) 
for x > 1 
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The limiting forms of scattering times in the small and 
large r s regimes are given by 

10c _ j [7r/2 + 12r s +8r s ln(r s ) + ...] for r s < 1 . - 
Ttc \ tt/32 - l/(15r.) + 0{l/r 2 ) for r s » 1 ( -> 

and 

roc _ ( r s /2 - nr 2 /2 - 2rf ln(r s ) + 
r sc \ 7r/16-l/(12r s ) + 0(l/r 2 ) 



for r s < 1 , . 
for r s > 1 v ; 



From Eq. (|8|) we can find easily the ratio of trans- 
port scattering time to the single particle relaxation time, 
Ttc/Tsc- Fig. Q] shows the ratio as a function of r s . Note 
that the ratio is only a function of r s and does not de- 
pend on the carrier density. We have the limiting form 
of Tt c /r sc , in the small and large r s regimes 
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- ^ln(r s ) for r s < 1 
for r s > 1 

57r r s 6 



(16) 



For large r s , the ratio approaches 2, which indicates that 
scattering is not isotropic even though the screening is 
strong. On the other hand, r tc /r sc diverges as l/r s for 
small r s . In r s — ► limits, where the screening is very 
weak, the transport scattering rate is much weaker than 
the single particle scattering rate, which indicates that 
the mobility in these limits is very high. 

For short-ranged ((5-function) scatterers with impurity 
density tin and potential strength vq we have the scat- 
tering times (considering screening effects) 
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In Fig. [T] we show the ratio of the scattering times for 
the short-ranged scatterers, T t s/r s s, as a function of r s , 
and the limiting forms are given by 



TtS 
TsS 



2 + 32r s ln(r s )/7r for r s < 1 
1 + 128r s /1057r for r s > 1 



(19) 
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We find that 1 < T tS /T sS < 2 (see Fig. [TJ. The ra- 
tio does not exceed 2 for short-ranged scatterers, but for 
charged impurity scattering the ratio is always greater 
than 2. Thus the ratio of the scattering times may of- 
fer the possibility of determining the relevant scattering 
mechanisms in disordered graphcnc layers. 

Now we compare the calculated scattering time ratio of 
graphene with the scattering time ratio of a normal 2D 
electron system with parabolic energy dispersion. We 
have the scattering time ratio for charged impurity (see 
Appendix A) in a normal 2D layer as 
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where r a = l/asy/im (&b is an effective Bohr radius of 
the system and n is a carrier density). Unlike graphene, 
in normal 2D systems r s depends on the carrier density. 
In small r s (or high density) limit the ratio shows the 
same leading order behavior as graphene (i.e. T t /r s oc 
l/r s ), but in large r s (low density) limit Tt/r s — > 1. The 
difference between normal 2D system and graphene in the 
strong screening limit can be traced back to the suppres- 
sion of the backscattering in graphene. We also find the 
difference in the scattering time ratio for short-ranged 
impurity. The ratio becomes 2/3 < T t s/T s $ < 1. (See 
Fig. [1] and Appendix A.) Thus, the ratio is always less 
than 1. In experiment, the ratio of the transport time to 
the single particle time has been found to be smaller than 
1 in Si-MOSFET systems. Thus, in Si-MOSFET system 
the short-ranged scattering (such as interface roughness 
scattering) dominates. 

It is, in principle, possible for the charged impurities to 
be at a distance u Zi" away from the 2D graphene layer. 
In fact, in 2D modulation doped GaAs-AlGaAs semicon- 
ductor heterostructures the charged dopants are put a 
distance d(= Zj) away (inside the GaAlAs insulating bar- 
rier region) from the 2D electron gas in order to minimize 
the degradation of the electron mobility due to remote 
dopant scattering. As discussed in Sec. I of this paper, 
a large separation (kpZi > 1) leads to a very strong en- 
hancement of the Tf/r s in the GaAs-GaAlAs heterostruc- 
ture since large-angle scattering by the remote impurities 
is strongly suppressed by the separation. To see the effect 
of separating the impurities from the graphene layer, all 
we need to do is to modify the form of the q-space Fourier 
transform of the electron-impurity Coulomb interaction 
by the factor e~ qZi arising from the separation between 
the 2D carriers in the graphene layer and the charged im- 
purity centers, leading to V e i(q, z,-) = V e i(q)e~ qZi . With 
this simple modification for V e i we can calculate the re- 
sults for Tt and r s for remote scatterers using the same 
formalism as above. 

In Figs. 1-4 we present our theoretical results for r t / t s 
as well as t 4 and t s individually in graphene comparing 
it with the corresponding regular 2D results (i.e. non- 
chiral 2D electron system with parabolic band disper- 
sion). We use exactly the same materials parameters 
for both graphene and the regular parabolic 2D system. 
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FIG. 1: (Color online) Calculated the ratio of the transport 
scattering time to single particle scattering time r t /r s (a) as 
a function of the parameter r s = e 2 /jk and (b) as a func- 
tion of k]p J qt f • Top (bottom) solid line represents the ratio 
for charged impurity (short-ranged neutral impurity) scatter- 
ing. Top (bottom) dashed line represents the ratio of nor- 
mal 2D system with parabolic dispersion for charged impurity 
(short-ranged neutral impurity) scattering. We note that for 
graphene kp/qTF = l/4r s . 



The results are plotted as a function of the dimension- 
less parameter r s (or kp/qTF), the impurity location 
from the interface z%, or the density n where kp (qtf) 
are respectively the Fermi wave vector (Thomas-Fermi 
wave vector) for the system. We give the expressions 
for Uf and qxF (with g as the spin and valley degen- 
eracy factor): kp = (Ann/g) 1 / 2 for both graphene and 
parabolic 2D systems, qTF — ge Icf /(kvf) for graphene 
and qxF = gme 2 / (kH 2 ) for 2D parabolic systems. We 
note that <?tf/^f = r s for g = 1 and ?tf/^f oc r s for 
any value of g for both graphene and parabolic 2D sys- 
tems. (In general, g — 4 for graphene.) 

The most important qualitative feature of Figs. 1-2 
is that in general rt/r s is larger in graphene than in the 
corresponding 2D parabolic system for all values of Zi in- 
cluding Zi = (i.e. when the impurities are right at the 
substrate-2D layer interface). For z^O and for r s w 1, 
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FIG. 2: (Color online) (a) The ration Tt/r a of graphene (with 
substrate Si02, corresponding to r s — 0.9) as a function of 
kF/qTF for different values of the separation Zj between the 
electron layer and the impurity layer. We use the dimen- 
sionless quantity ZiqrF = 0, 0.1, 0.5. Solid (dashed) lines 
indicate results for graphene (normal 2D system), (b) The 
same as Fig. 2(a) as a function of z% for two different electron 
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FIG. 3: (Color online) (a) Scattering times (in ps) of 
graphene (with substrate Si02) as a function of electron den- 
sity for two different impurity locations 2, = 0, 10A. r tc (r sc ) 
indicates the transport scattering time (single particle relax- 
ation time) for Coulomb impurity scattering. We use an im- 
purity density m c = 10 12 cm~ 2 . The same as Fig. 3(a) for 
a neutral impurity scattering. We use an impurity density 



1, 5 x 10 cm" 2 . Solid (dashed) line denote the n iS = 2 x 10 cm and v Q = IKeVA 2 



densities, 

results for charged Coulomb impurity (short-ranged neutral 
impurity). 



the situation of current interest in graphene where the 
impurities are thought to be within ~ 1 nm of the inter- 
face (and the substrate is usually SiC^), we get r t /T s ~ 2 
in graphene although it is expected to increase fast with 
increasing Zi. We therefore suggest a measurement of 
r t /t s as a spectroscopic tool for the determination of the 
location of the impurity centers in graphene. It is inter- 
esting to note that in general {T t /T s ) graphene > (T t /T S ) 2 D 
with the difference between them increasing sharply with 
increasing kp/qTF or Zj. 

In Fig. 3 we show our calculated graphene transport 
and single particle times individually as a function of car- 
rier density for various values of the impurity location z^. 
In Fig. 3(a) where charged impurity Coulomb scattering 
is involved, both Tt and r s increase with carrier density, 
but Tt increases much faster, making T t /r s increase with 



increasing density. But for short-range scattering (Fig. 
3(b)), the reverse is true with both Tt and r s decreasing 
with increasing density, but t s decreasing faster, again 
leading to an increasing Tt /t s with increasing carrier den- 
sity. 

Finally, in Fig. 4 we show our calculated level broad- 
ening r = h/2r s as a function of carrier density in 
graphene. We note that for charged impurity scatter- 
ing T/Ep decreases monotonically with increasing carrier 
density. However, the scaled damping rate due to short- 
ranged impurity scattering is independent of the carrier 
density. 



III. DISCUSSION AND CONCLUSION 

We calculate the transport scattering time (r t ) and 
the single particle relaxation time (r s ) for disordered 
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FIG. 4: (Color online) Calculated damping rates scaled by 
Fermi energy Y / Ef as a function of carrier density for r s — 0.9 
which corresponds to graphene sample with Si02 substrate. 
Solid lines indicate damping rates due to charged impurity 
scattering with impurity density rii C = 100, 10, 1 x 10 10 cm~ 2 
(from top to bottom), respectively. Dashed line indicate a 
damping rate due to short-ranged impurity with impurity 
density na = 10 10 cm -2 and potential strength vo —1 KeV 
A 2 . We define Y = h/2r s . 



APPENDIX A 

Here we provide the scattering times of a normal 2D 
system with parabolic band. Using 2D RPA screening 
function^ at T = we have the scattering times for 
charged impurity centers 



1 

r tc 
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where 
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(Ala) 
(Alb) 
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rii C is the density of charged impurity, qo — qTF/^kp 
(qtf — g/o-B is a 2D Thomas- Fermi wave vector with 
effective Bohr radius <xb = h 2 /me 2 ), and Itc(qo)> -TscGzo) 
are given by 



Itdqo) 



(A3) 



graphene in the lowest order of the electron-impurity in- 
teraction (Born approximation). The scattering mecha- 
nisms which we consider are screened charged impurity 
scattering and short-range potential (e.g. caused by lat- 
tice defects). For the screening function we use the RPA. 
We find that for short ranged scatterers, the ratio of the 
scattering time to the single particle relaxation time is 
always smaller than 2, but for charged impurity scatter- 
ing, the ratio is always greater than 2. These theoretical 
results provide a technique to experimentally discrimi- 
nate between short-range and Coulomb scattering as to 
the relevant scattering mechanism in disordered graphene 
layers. We also find a strong dependence of the ratio r t /r s 
in graphene on Zi the separation of the impurities from 
the 2D graphene layer. Somewhat surprisingly the de- 
pendence of Tt /t s on Zi is stranger in graphene than in the 
corresponding parabolic 2D system, leading to the possi- 
bility that an accurate measurement of Tt /r s in graphene 
could lead to better understanding of the impurity loca- 
tion underlying graphene disorder. In particular, T t /r s 
increases rapidly with increasing Zj in graphene (similar 
to non-chiral 2D GaAs-GaAlAs modulation doped het- 
erostructure) and this could be directly experimentally 
tested. 

We conclude by emphasizing that independent mea- 
surements of Tt and t s in graphene samples could lead to 
detailed useful insight into the nature of disorder scatter- 
ing of graphene carriers. 

Acknowledgments 

This work is supported by U.S. ONR, NSF-NRI, and 
SWAN. 



Isc{qo) = -d^ fiqQ) 



(A4) 



where f(x) is given in Eq. (|12[) and its derivative is given 
by 



df 1 1 
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For short-ranged impurity centers we have 
1 1 
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(A5) 

(A6a) 
(A6b) 

(A7) 



riis is the density of short-ranged impurity, m$ , vq is the 
potential strength and Its(qa), I s $(qo) are given by 

Its(qo) = J - 4<? + 'Wo - 2/- [q 4 f(qo)] (A8) 
Z aqo 



Isc{qo) = 7; - -r- [QofiOo)] 
Z aq 



(A9) 
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